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Abstract. Although the first millisecond pulsars (MSPs) were discovered 30 years ago we still 
(N. do not understand all details of their formation process. Here, we present new results from 
Tauris, Langer & Kramer (2012) on the recycling scenario leading to radio MSPs with helium 
or carbon-oxygen white dwarf companions via evolution of low- and intermediate mass X-ray 
binaries (LMXBs, IMXBs). We discuss the location of the spin-up line in the PP-diagram and 
estimate the amount of accreted mass needed to obtain a given spin period and compare with 
observations. Finally, we constrain the true ages of observed recycled pulsars via calculated 
isochrones in the PP-diagram. 
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1. Introduction 

Q-i Binary MSPs represent the advanced phase of stellar evolution in close, interact- 

Q ■ ing binaries. Their observed orbital and stellar properties are fossil records of their 
I evolutionary history. Thus one can use binary pulsar systems as key probes of stel- 
C/3 . lar astrophysics. Although the standard recycling scenario (Alpar et al. 1982; Bhat- 
tacharya & van den Heuvel 1991) is commonly accepted, many aspects of the mass- 
transfer process and the accretion physics are still not understood in detail. Examples of 
such ambiguities include the accretion disk structure, the disk-magnetosphere transition 
zone, the accretion efRciency, the decay of the surface B-field of the neutron star and the 
outcome of common envelope evolution. For further details on these aspects, details in 
■ general and discussions of our results, we refer to our journal paper, Tauris et al. (2012). 

o; 

2. The spin-up line 

• Some of the above mentioned simplifications become a problem virhen trying to probe 

the formation and the evolution of observed recycled radio pulsars located near the 
classical spin-up line for Eddington accretion in the PP-diagram (e.g. as illustrated with 
the MSP J1823-3021A, Freire et al. 2011). The location of the spin-up Hue can be found 
^ . by considering the equilibrium configuration when the angular velocity of the neutron 
- - ' star is equal to the Keplerian angular velocity of matter at the magnetospheric boundary 
(r,„ag) where the accreted matter enters the magnetosphere, i.e. il^, = fl^q — lOq f^K(^mag) 
or: Poq = 27r(rj^jjjg/GAf)^/^ oj"^, where 0.25 < ^ 1 is the so-called critical fastness 
parameter. Introducing the magnetospheric coupling radius, (j) = r-mag/'^Aifvcn and the 
magnetic inclination angle, a we can rewrite this expression: 

P = + '''' ^^-^^ 
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which can be plotted directly in the PP-diagram. (Here M is the mass of the pulsar, M 
is its accretion rate and / is its moment of inertia.) In case sin a = ^ = Wc = 1 we find: 

P = 3.7 X 10-19 (M/Mo)V3 P^/3 (JL.) (2.2) 

where P^g is the equilibrium spin period in units of milliseconds. We have included the 
plasma term in the spin-down torque using the combined model of Spitkovsky (2006) to 
compensate for the incompleteness of the vacuum magnetic dipole model. 

In Fig. [T] we have plotted equation (12. ip for different values of a, (f) and to illustrate 
the uncertainties in the applied accretion physics to locate the spin-up line. In all cases 
we assumed a fixed accretion rate of M = M-^dd- The location of the spin- up line is 
simply shifted one order of magnitude in P down (up) for every order of magnitude M 
is decreased (increased). It is important to realize that there is no universal spin- up line 
in the PP-diagram. Not only are M and a individual to each pulsar (giving rise to the 
width of each band), also and uic could be related to B, a and M. 
If we assume that accretion onto the neutron star is indeed Eddington limited, then the 
three bands in Fig. [T] represent upper limits for the spin-up line for the given sets of 
((f), ujc)- Thus we can in principle use this plot to constrain ((/>, ujc) and hence the physics 
of disk-magnetosphere interactions from future detections of MSPs. 




Spin period, P (ms) 

Figure 1. Calculations of three spin-up lines, shown as coloured bands, depending on the 
parameters (0, ljc)- The upper boundary of each band (or "line") is calculated for a neutron 
star mass M = 2.0 M© and a = 90°. The lower boundary is calculated for M = 1.0 M© and 
a — 0° . The green (central) hatched band corresponds to (j) = 1 and ojc = 1. The blue and 
red hatched bands are upper and lower limits set by reasonable choices of the two parameters 
{(p, oJc)- In all three cases the spin-up line is calculated assuming accretion at the Eddington 
limit, M = MEdd- The observed distribution of binary and isolated radio pulsars in the Galactic 
disk are plotted as filled and open circles, respectively. Also plotted is the pulsar J1823— 3021A, 
located in the globular cluster NGC 6624. (Fig. adapted from Tauris et al. 2012.) 
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When modeling the birth spins of recycled radio pulsars it is important to include the 
braking torque acting during the Roche- lobe decoupling phase (RLDP) when the donor 
star terminates its mass transfer. It has been shown that accreting X-ray MSPs may lose 
up to 50 % of their rotational energy during the RLDP of LMXBs (Tauris 2012). 

3. Relation between accreted mass and final spin period 

Recycled pulsars obtain their fast spins via angular momentum exchange from the 
differential rotation between the accretion disk and the neutron star. The amount of spin 
angular momentum added to an accreting pulsar is given by: 

A J, = j n{Lu, t) M{t) ^GM{t)r^,^{t)i{t) dt (3.1) 

where n{uj,t) is a dimensionless torque. Assuming n{uj,t) = 1, and M{t), r,„ag(i) and 
^ (t) to be roughly constant during the major part of the spin- up phase we can obtain a 
simple and convenient expression to relate the (minimum) amount of accreted mass and 
final equilibrium spin period (see also Alpar et al. 1982): 

AM., = 0.22M,M^ ,3.2) 

Puis 

assuming a numerical factor f{a,^,(f>,ujc) — 1 (from disk-magnetosphere interactions). 
Considering a pulsar with a final mass of IAMq and a recycled spin period of either 
2 ms, 5 ms, 10 ms or 50 ms requires accretion of 0.10 M©, 0.03 Mq, 0.01 Mq or lO^^Mg, 
respectively. Therefore, it is no surprise that observed recycled pulsars with massive 
companions (CO/ONeMg WD or NS) in general are much more slow rotators - compared 
to MSPs with He WD companions - since the progenitor of their massive companions 
evolved on a relatively short timescale in IMXBs or HMXBs, only allowing for very little 
mass to be accreted by the pulsar. 

4. True age isochrones of recycled pulsars 

In order to investigate if we can understand the distribution of MSPs in the PP- 
diagram we have traced the evolution of eight hypothetical, recycled MSPs with different 
birth locations {Pq, Pq). In each case we traced the evolution as a function of age, t for 
a constant braking index 2 ^ n ^ 5 and calculated isochrones by integration for each 
pulsar given that P{t, n, Pq, Pq). The results are shown in Fig. [2] together with observed 
data. Furthermore, we plotted two isochrones (see fat green and pink lines) calculated 
for {Pq — 1.0 ms, n — 3, t — 1.5 Gyr) and {Pq — 7.0 ms, n = 3, i = 12Gyr), respectively, 
and with no restrictions on Pq (or Bq). 

A number of interesting conclusions can be drawn from this diagram. The overall distri- 
bution of observed pulsars follows nicely the banana-like shape of the two fat isochrones, 
see also Kiziltan & Thorsett (2010), and hence pulsars are recycled with a wide range 
of final B-ficlds. Although these curves are not an attempt for a best fit to the obser- 
vations it is interesting to notice that close to 90 % of all recycled pulsars (even up to 
P = 100 ms) are compatible with being born (recycled) with an initial spin period of 
Pq = 1 — 7 ms and having ages between 1.5 and 12 Gyr. However, from a binary evolution 
point of view many of the P = 20 — 100 ms pulsars (the mildly recycled pulsars) are 
born with such relatively slow spins (see Section and hence they need not be that old. 
This can, for example, be verified by cooling age determinations of their WD companion 
stars. Pulsars with small values of the period derivative, P c± 10^^^ hardly evolve at all in 
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Figure 2. Isochrones of eight iiypotlietical recycled pulsars born at the locations of the red 
stars. The isochrones were calculated for different values of the braking index, 2 ^ n 5C 5. Also 
plotted are inferred B- field values (dashed lines) and characteristic ages, r (dotted lines). The 
thin gray lines are spin-up lines with Af/MEdd = 1, 10~^, 10~^, 10"'^ and 10~* (top to bottom, 
and assuming sin a = <^ = cjc = !)• In all calculations we assumed a pulsar mass of 1.4 M©. It 
is seen how the banana shape of the two fat isochrones (see text) fits very well with the overall 
distribution of observed pulsars in the Galactic disk. Binary pulsars are marked with solid circles 
and isolated pulsars are marked with open circles, using data from the ATNF Pulsar Catalogue 
and corrected for kinematic (Shklovskii) effects. (Fig. adapted from Tauris et al. 2012.) 

the diagram over a Hubble time. This trivial fact is important since it tells us that these 
pulsars were basically born with their currently observed values of P and P (first pointed 
out by Camilo et al. 1994). Hence, some pulsars with characteristic ages of r ~ 100 Gyr 
could in principle have been recycled very recently ~ demonstrating the unreliability of 
r as a true age indicator. It is also interesting to notice PSR J1801— 3210 (discovered 
by Bates et al. 2011) which must have been recycled with a relatively slow birth period, 
Po ~7 ms despite its low B-field < 10^ G, see Fig. [1] for its location in the PP-diagram. 
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